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Abstract 

We study magnetically charged black holes in the Einstein- Yang-Mills-Higgs theory in the 
limit of infinitely strong coupling of the Higgs field. Using mixed analytical and numerical meth- 
ods we give a complete description of static spherically symmetric black hole solutions, both 
abelian and nonabelian. In particular, we find a new class of extremal nonabelian solutions. 
We show that all nonabelian solutions are stable against linear radial perturbations. The im- 
plications of our results for the semiclassical evolution of magnetically charged black holes are 
discussed. 
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1 Introduction 

It has long been thought that the only static black hole solution in spontaneously broken 
gauge theories coupled to gravity (Einstein- Yang-Mills-Higgs (EYMH) theories) is the Reissner- 
Nordstr0m (RN) solution, with covariantly constant Higgs field and the electromagnetic field 
trivially embedded in a nonabelian gauge group [1]. This belief, based on no-hair properties of 
black holes, was put in doubt by the discovery of essentially nonabelian black holes in EYM 
theory (so called colored black holes) [2], because it became clear that the uniqueness properties 
of black holes in Einstein-Maxwell theory [3] are lost when the gauge field is nonabelian. This 
suggested that there might also exist nontrivial (i.e. different than RN) black holes when the 
Higgs field is included. Indeed, it was found recently that in the EYMH theory there are spher- 
ically symmetric magnetically charged black holes which are asymptotically indistinguishable 
from the RN solution but carry nontrivial YM and Higgs fields outside the horizon [4,5]. These 
essentially nonabelian solutions may be viewed as black holes inside 't Hooft-Polyakov magnetic 
monopoles. Because the EYMH field equations are rather complicated (even for spherical sym- 
metry), one does not have yet a complete picture of these solutions. In particular, one would like 
to know how the solutions depend on external parameters (like coupling constants)^], whether 
there are solutions with degenerate horizon, and, most important, whether the solutions are 
stable. 

The aim of this paper is to give answers to these questions within a model obtained from the 
EYMH system by taking the limit of infinitely strong coupling of the Higgs field. In this limit the 
Higgs field is effectively "frozen" in its vacuum expectation value. Thus one degree of freedom 
is reduced which simplifies considerably the analysis of solutions. This model is interesting on 
its own for testing ideas about black holes and as we will argue has the additional virtue that 
many of its features carry over to the full EYMH theory (albeit there are also some important 
differences) . 

The paper is organized as follows. In the next section we define the model and derive the 
spherically symmetric field equations. 

In Section 3 we consider static black hole solutions. The analysis of asymptotic behaviour 
of solutions at infinity and at the horizon allows us to obtain information about possible global 
1 This question was examined in the Prasad-Sommerfield limit in Ref. [5]. 
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behaviour of solutions. We show that black hole solutions may exist only if their parameters 
satisfy certain inequalities. In particular, we derive necessary conditions for the existence of 
extremal black holes. 

In Section 4 we describe numerical results and discuss the qualitative properties of solutions. 
The bifurcation structure of solutions as functions of external parameters is discussed in detail. 
We also comment on the status of no-hair conjecture in our model. 

Section 5 is devoted to stability analysis. We derive the pulsation equation governing the 
evolution of small radial perturbations about a static solution and solve it numerically. Special 
emphasis is placed on the relation between stability properties and the bifurcation structure of 
solutions. Using this relation we conjecture about the stability of black holes in the full EYMH 
theory. 

Finally, in Section 6 we summarize our results and touch upon the question of evolution of 
black holes in the model due to the Hawking radiation. 



2 Field Equations 

The Einstein- Yang-Mills-Higgs (EYMH) theory is defined by the action 
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F = dA + e[A,A] 
Dcj) = d<p + e[A,4>] 



(1) 

(2) 
(3) 
(4) 



where the YM connection A and the Higgs field (j> take values in the Lie algebra of the gauge 
group Q. Here we consider Q = SU{2). 

There are four dimensional parameters in the theory: Newton's constant G, the gauge cou- 
pling constant e, the Higgs coupling constant A and the vacuum expectation value of the Higgs 
field v. Out of these parameters one can form two scales of length: 

length. 
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The ratio of these two scales 

a = VGv (5) 
plays a role of the effective coupling constant in the model. The second dimensionless parameter 

measures the strength of the Higgs field coupling. In this paper we consider in detail the case 
f3 = co. In this limit the Higgs field is "frozen" in its vacuum expectation value v. 

We are interested in spherically symmetric configurations. It is convenient to parametrize 
the metric in the following way 

ds 2 = -B 2 Ndt 2 + N~ l dr 2 + r 2 ((M 2 + sin 2 tfdip 2 ) (7) 

where B and N are functions of (t,r). 

We assume that the electric part of the YM field vanishes. Then the purely magnetic spher- 
ically symmetric su(2) YM connection can be written, in the abelian gauge, as 

eA = WTid'd + (cot #T3 + WT2) sin fldip (8) 

where Tj (i = 1, 2, 3) are Pauli matrices and w is a function of (t, r). 

The spherically symmetric Ansatz for the Higgs field in the abelian gauge is 

(j) = VT 3 h(t,r). (9) 

As we said we concentrate on the limit (3 = 00, which implies h(r,t) = 1. Note that although 
the Higgs field is constant, its covariant derivative is nonzero and generates a mass term for the 
YM field: 

(Dcj)) 2 = v 2 w 2 . 

The YM curvature is given by 

eF = wdt Afi + w'dr A n - (1 - u> 2 )r 3 cM A sin tidip (10) 
where (' = dt, ' = d r ) and 

$7 = Tidfl + T2 sin -ddip. 
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From the EYMH Lagrangian (1) one derives the following components (in orthonormal frame) 
of the stress energy tensor T a ^ 

/> = i T « = ^ Nw ' 2 + b ~ 2n ~^ + I w ! + " 2 S <"> 



"=SF T ^ = 2fl "'i W - (13> 



(In what follows we do not make use of the other nonzero components = T^.) 
The Einstein equations reduce to the system 

N' = -(1 -N)- 2Grp (14) 

B' = GrBN- 1 {p-r) (15) 
TV = -2GrBNa. (16) 
The remaining #$-component of the Einstein equations is equivalent to the YM equation 

- (B^N^wJ + (BNw'Y + ^Bw(l - w 2 ) - Be 2 v 2 w = 0. (17) 

For G = the equations (14-16) are solved either by the Minkowski spacetime N = B = 1, or 
by the Schwarzschild spacetime = 1, N = 1 — 2m/r, and then the Eq. (17) reduces to the 
YMH equation on the given background. For v = 0, the system (14 - 17) reduces to the EYM 
equations. 

3 Static Solutions 

In this section we wish to consider static black hole solutions of the system (14 - 17). In terms 
of the dimensionless variable x = the static equations are 

N> = 1(1 - N) - 2 -(N W ' 2 + + a 2 W 2 ) (18) 

B' = -Bw' 2 (19) 



(NBw'Y + b\w{1 - w 2 - a 2 x 2 ) = 0. (20) 
ar 
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Note that the function B can be eliminated from Eq. (20) by using Eq. (19). This system of 
equations has been studied previously by Breitenlohner et al. [5] who found asymptotically flat 
solutions with naked singularity at x = 0. This singularity is an artefact of the limit (3 = oo, 
since then there is no symmetry restoration at the origin and the last term in the expression 
(11) for the energy density diverges at r = 0. In the black hole case the singularity is hidden 
inside the horizon. 

We will consider solutions of Eqs. (18 - 20) in the region x 6 [x#,oo), where xu is the 
radius of the (outermost) horizon. The boundary conditions at x = xu are 

N(x H ) = 0, N'(x H ) > 0, (21) 

and the functions w, w' , B are assumed to be finite. 

At infinity we impose asymptotic flatness conditions which are ensured by 

iV(oo) = 1, w{oo) = w'(oo) = 0, B{oo) finite. (22) 

One explicit solution satisfying these boundary conditions is well known. This is the Reissner- 
Nordstr0m (RN) solution 

2m 1 

w = 0, B = 1, N = 1 + (23) 

x x z 

where mass m > 1 (mass is measured in units 1/y/Ge). The corresponding YM curvature is 

F = —Tsdi} A smfidip 

showing that it is an abelian solution describing a black hole with unit magnetic charge (the 
unit of charge is 1/e). 

It turns out that the system (14-16) admits also nonabelian solutions which we have found 
numerically. Before discussing the numerical results we want to make some elementary observa- 
tions about the global behaviour of solutions satisfying the boundary conditions (21) and (22). 
In what follows we assume that a is nonzero. 

First, consider the behaviour of solutions at infinity. The asymptotic solution of Eqs. (18 - 
20) is 

w ~ e~ ax 

N ~ 1 - — + \ + 0(e- 2ax ) 
x x z 

B ~ l + 0{e~ 2ax ) 
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hence for large x the solutions are very well approximated by the RN solution (23). 

Second, notice that w has no maxima for w > 1 and no minima for w < — 1, which follows 
immediately from Eq. (20). Thus, if w once leaves the region w G (—1,1), it cannot reenter 
it. Actually, w stays within this region for all x > xh, because w(xh) < 1 (without loss of 
generality we may assume that w(xh) > 0, because there is a reflection invariance w — > —w). 
To see this, consider Eq. (20) at x = xh- Assuming that w"(xh) is finite one has 



N'w' + \w{\ - w 2 - a 2 x 2 ) 
x z 



= 0. (24) 

X=X H 



If w(xh) > 1, then w'(xh) > 0, but since there are no maxima for w > 1 the condition w(oo) = 
cannot be fulfilled. Thus w(xh) < 1- Next one can show that w has no positive minima and 
negative maxima. Suppose that there is a positive minimum at some xq > xh- Then Eq. (20) 
implies that the function f(x) = 1 — w 2 — a 2 x 2 is negative at xq. For x > xq, f(x) decreases 
(because w' > 0), hence there cannot exist a maximum of w for x > xq, so again the condition 
w(oo) = cannot be met. One can repeat this argument to show that w'(xh) < 0, because 
if w'(xh) > 0, then /(xh) < 0, as follows from Eq. (20). Finally, notice that for x > l/a, 
w cannot have any extrema because then f(x) < 0, hence there are no positive maxima and 
negative minima (whereas other extrema we have already excluded). 

To summarize, we have shown that if there exits a solution of Eqs. (18 - 20), satisfying the 
boundary conditions (21-22), then the function w stays in the region w G (—1,1) and either 
monotonically tends to zero or oscillates around w = for x < l/a and then monotonically 
goes to zero. 

Now, we will show that black hole solutions may exist only if the parameters a and xh 
satisfy certain inequalities. Let b = w(xh)- We have shown above that w'(xh) < for b > 0, 
hence from (24) we obtain 

1 - b 2 

-^T- > « 2 ( 25 ) 
x H 

which implies the necessary condition for the existence of a black hole solution 

ax H < 1. (26) 

As we shall see in the next section this is not a sufficient condition. 

For xh < 1 we can improve the inequality (26). From Eq. (18) we have 

N . iXH) = ±-^-^ (27) 

XH X H XH 
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hence 

(1 - b 2 ) 2 -x 2 H + 2a 2 x 2 H b 2 < 0. (28) 
This inequality has a real solution for b only if the discriminant 

U = {a'x 2 H - 2a 2 + l)x 2 H (29) 

is nonnegative, which is always satisfied for xh > 1 but for xh < 1 this implies the additional 
condition 

a 2 <a 2 max = ±-(l-yjl-x%). (30) 
x H 

It turns out from numerical results that the inequality (30) is also a sufficient condition for the 
existence of black holes. Note that when xh — ► then a max — > 1/ V2, which is the upper bound 
for "regular" solutions found by Breitenlohner et al. [5]. 

Finally, let us see whether there may exist extremal black holes in the model. By extremal 
we mean a solution with degenerate horizon, i.e. N'(xh) = 0. For such solutions the inequalities 
(25) and (28) are saturated. Eliminating b 2 we obtain the condition 

a^x\ - 2a 2 + 1 = (31) 

which is equivalent to A = and is solved by a = a max given by (30) provided that xh < 1- Let 
us point out that the necessary conditions for the existence of extremal black hole solution in 
EYMH theories were derived some time ago by Hajicek [6]. In the terminology of Hajicek, Eq. 
(31) is a special case of the zeroth-order condition. In the next section we will find numerically 
the extremal solutions satisfying (31). 



4 Numerical Results 



The formal power-series expansion of a solution near the horizon is 

w(x) = b+J2yW ik \x H )(x-x H )\ 

k=l 

oo -. 

k=l K - 

All coefficients in the above series are determined, through recurrence relations, by 6, in particu- 
lar the expressions for w'(xh) and N'(xh) are given by Eqs. (24) and (27). Thus this expansion 
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(assuming that its radius of convergence is nonzero) defines a one-parameter family of local 
solutions labelled by the initial value b. We use a standard numerical procedure, called shoot- 
ing method, to find such values of b for which the local solution extends to a global solution 
satisfying the asymptotic boundary conditions (22). Note that b is not arbitrary but, as follows 
from Eqs. (25) and (28), must lie in the interval 

1 - o?x\ - x H \loi A x 2 H - 2q 2 + 1 < b 2 < 1 - o?x\. (33) 

We find that for every xh , there is a maximal value o; ma x (xh)i such that for a > owx(^tf ) there 
are no solutions, while for ex ^ 0! max there is exactly one solution for which w is monotonically 
decreasing. This solution we call fundamental, in contrast to solutions with oscillating w which 
we call excitations. In what follows we restrict our attention to fundamental solutions. At the 
end we will briefly describe excitations. 

The numerical results fo several values of xh are summarized in Table 1. When xh < 1, 
then a max (x#) is given by Eq. (30). For xh > 1, amax(sff) is displayed in Table 2 (in this case 
the analytical bound given by Eq. (26) is not sharp). In Fig. 1 we plot w{x) for xh = 2 and 
several values of a. 

Table 1: Shooting parameter b and mass m (in units \ j\fGe) as functions of a for xh = 0.5, 
xh = 1, and xh = 2. 



xh = 0.5 


xh = 1 


x H = 2 


a 


b 


m 


a 


b 


m 


a 


b 


m 


0.01 


0.99993 


0.268 


0.01 


0.99975 


0.517 


0.01 


0.99907 


1.018 


0.2 


0.98296 


0.572 


0.1 


0.98191 


0.663 


0.1 


0.9317 


1.147 


0.4 


0.94782 


0.809 


0.5 


0.63099 


0.9805 


0.2 


0.7142 


1.2286 


0.6 


0.923028 


0.954 


0.74 


0.1852 


0.99937 


0.26 


0.43074 


1.248 


0.7 


0.928359 


0.990 


0.9 


0.00868 


0.99999 


0.28 


0.2414 


1.2498 


0.73 


0.930534 


0.996 








0.28867 


0.012 


1.25 
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Table 2: Maximal value of a as a function of xh for xh > 1. 



XH 




1 


1 


1.000001 


0.822 


1.01 


0.744 


1.1 


0.605 


1.5 


0.398 


2 


0.288 


10 


0.055 


20 


0.027 



On the basis of analytical and numerical results we have a pretty clear picture of the qual- 
itative behaviour of solutions. Let us discuss in more detail how the solutions depend on the 
parameters a and xh- First consider the limit a — > 0. In our units this corresponds to v — > 0. 
In this limit a solution on any finite region outside the horizon is well approximated by the 
Schwarzschild solution 

1 AT 1 XH 

w = 1, N = 1 , 

with mass m = xjj/2. However, for large x and small but nonzero a the term a 2 x 2 in Eqs. (18) 
and (20), becomes dominant and asymptotically the solution tends to the RN solution. 

Next, consider the limit a — > a max (xij). In this case the behaviour of solutions depends 
on whether xh is less or greater than one. For xh < 1, the interval of allowed values of 
b, given by Eq. (33), shrinks to zero as a goes to a max , and for a = a max , b is uniquely 
determined by xh- The corresponding limiting solution describes an extremal black hole which 
is essentially nonabelian. For xh > 1, when a goes to a max the solution tends to the RN 
solution and for a = a max coalesce with it. Thus for given xh > 1, the point a = a max (a;#) 
is a bifurcation point: for a > a max there is only one solution (RN), while for a < a max the 
second (nonabelian) solution appears. The bifurcation diagram in the plane (a,xn) is graphed 
in Fig. 2. The nonabelian solutions exist only in the region below the curve ABC. Along the 
curve AB, the nonabelian solutions are extremal. Along the curve BC the nonabelian and RN 
solutions coalesce. The RN solutions exist for > 1 and do not depend on a (for xh = 1 the 
RN solution is extremal). 
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The bifurcation of solutions is also shown in Fig. 3, where the mass m as a function of a is 
graphed for given xh > 1- Notice that the mass of the RN solution ttirn = \ (%H + ^) is larger 
than the mass of the nonabelian solution with the same xh- This suggests that for a < a max , 
where there are two distinct solutions with the same radius of the horizon xh, the RN solution 
is unstable. We will show in the next section that this is indeed the case. 

Let us comment on the status of the no-hair conjecture in our model. The no-hair conjecture 
(in its strong version) states that stationary black hole solutions, within a given model, are 
uniquely determined by global charges defined as surface integrals at spatial infinity such as 
mass, angular momentum and electric or magnetic charge. In our case, all solutions are static 
and have unit magnetic charge, so the only global parameter by which the solutions may differ 
at infinity is their mass. For a > 1 there is only one solution for given mass m, namely the 
RN solution, hence the strong no-hair conjecture is valid. However, for a < 1 the situation 
is different. This is illustrated in Table 3 which shows how the masses of nonabelian and RN 
solutions depend on xh for given a. 

Table 3: Masses of nonabelian and RN black holes as functions of xu for a = 0.1 



X H 




"1-RN 


0.1 


0.226 




0.2 


0.274 




1 


0.663 


0.75 


2 


1.147 


1.25 


5 


2.599 


2.6 


5.55 


2.86508 


2.86509 


5.595 


2.886865 


2.886865 


6 




3.08333 


8 




4.0625 



For the RN solution xh > 1, so its mass ttt-rn = 9 {xh + ^— ) is bounded from below by 0.75. 
At x H &x {a) the nonabelian and RN solutions coalesce, so the maximal mass of the nonabelian 
solution is equal to m^(x H ax ). For given a < 1 and 0.75 < m < m^(x H ax ) there are two 



11 



distinct (i.e. RN and nonabelian) solutions with the same massQ. This violates the strong no-hair 
conjecture. 

We will show in the next section that in the region where two distinct solutions coexist, only 
one of them is stable. Thus, in our model, a weak no-hair conjecture holds, i.e. global charges 
determine uniquely the stable black hole solution. However, as we will argue, in the full EYMH 
theory even the weak no-hair conjecture is violated. 

Finally, let us briefly consider solutions for which w is not monotonic. Such solutions may 
be viewed as excitations of the fundamental solutions described above. Their existence in the 
EYMH model was first noticed in Ref. 5. To understand the existence of excitations it is useful 
to consider the limit a — > 0. Then, the Eqs. (18 - 20) become the EYM equations, which 
have a countable family of so-called colored black holes [2]. These solutions are labelled by 
an integer n = number of nodes of w, and in contrast to the a ^ case, they have zero 
magnetic charge (because w(oo) = ±1). For sufficiently small a, the excitations may be viewed 
as singular perturbations (in a) of colored black holes. As long as ax <C 1, the excitations 
are well approximated by colored black hole solutions (with the same number of nodes for w), 
but for large x the term a 2 x 2 in Eqs. (18) and (19) becomes dominant and forces w to decay 
exponentially. This is shown in Fig. 5, where the n = 1 excitation is graphed for two values of 
a and compared to the n = 1 colored black hole solution. 

It is rather difficult to find numerically the excitations with large n, however we expect that 
for given xh and the number of excitations is finite (in contrast to the a = case, where 

there are infinitely many solutions). Our expectation is based on the monotonicity of solutions 
for x < 1/a, proven in Section 3. The existence of excitations with arbitrarily large n would be 
inconsistent with this property, because, in analogy to the colored black holes, the location of 
nodes is expected to extend to infinity as n increases. 

The bifurcation structure of excitations is similar to that for the fundamental solution. For 
given xh > 1 there is a decreasing sequence {aSLtx; a max> • • •} °f bifurcation points such that 
at a^ax the n-th excitation bifurcates from the RN branch (by zeroth excitation we mean the 
fundamental solution). As we will discuss in the next section this picture is closely related to 
the stability properties of the RN solution. 



For sufficiently small a the excitations appear so there are even more than two solutions with the same mass. 
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5 Stability Analysis 

In this section we address the issue of linear stability of the black hole solutions described above. 
To that purpose we have to study the evolution of linear perturbations about the equilibrium 
configurations. Since we do not expect nonspherical instability (and since, admittedly, the anal- 
ysis of nonspherical perturbations would be extremely complicated), we restrict our study to 
radial perturbations. For radial perturbation the stability analysis is relatively simple, because 
the spherically symmetric gravitational field has no dynamical degrees of freedom and therefore 
the perturbations of metric coefficients are determined by the perturbations of matter fields. 
This was explicitly demonstrated for the EYM system by Straumann and Zhou [7] who derived 
the pulsation equations governing the evolution of radial normal modes of the YM field. Below 
we repeat their derivation with a slight modification due to the presence of the mass term in 
Eqs. (14 - 17). 

We define the functions a and b by e a = BN and e b = N and write the perturbed fields as 

w(r,t) = wo(r) + 5w(r,t) 

a(r,t) = a (r)+Sa(r,t) (34) 
b(r,t) = bo(r) + 6b(r,t) 

where (wo,ao,bo) is a static solution. We insert the expressions (34) into Eqs. (14 - 17) and 

keep first order terms in the perturbations. Hereafter we use dimensionless variables 

e . e 

t and x = — = r 



'G VG 
and we omit the subscript for static solutions. From (16) we obtain 



4 

5b = w'5w. (35) 

x 



The asymptotically flat solution of this equation is 



Linearization of Eq. (15) yields 
Thus, using (36), we obtain 



4 

5b = — w'5w. (36) 

x 



5a' - 5b' = -w'5w'. (37) 

x 



4 4 

5a' = -^w'5w - -w"5w. (38) 

ar x 
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Multiplying Eq. (17) by e a and linearizing we get 



- e~ 2a 5w + e- a (e a 5w'y + 6a' vf - 5be~ b w 



r (l-„,2)_ a S 



+ e 



1 



(1 - 3w 2 ) - a 2 



Now, we insert (36) and (38) into (39) and make the Ansatz 



to obtain the eigenmode equation 



S W = e^ T <£(x) 



where 



e~ a U = 5 e a (l + ax)w l 



x* 



-e a - b ww' 

x 



a 



„a— 6 



-o (1 - 3w 2 



a 



5w = 0. 
(39) 

(40) 
(41) 
• (42) 



The potential U is a smooth bounded function which vanishes at xh and tends to a 2 for x — > oo. 

One can introduce the tortoise radial coordinate to transform Eq. (41) into the one- 
dimensional Schrodinger equation. However we will not do so, because the numerical analysis 
is easier when one uses the coordinate x. 

A static solution (w, a, b) is stable if there are no integrable eigenmodes £ with negative to 2 . 
To check this we have applied a (slightly modified) rule of nodes for Sturm-Liouville systems 
[8], which states that the number of negative eigenmodes is equal to the number of nodes of the 
zero eigenvalue solution. Namely we have considered Eq. (41) with negative co 2 and looked how 
the solution satisfying £(xjj) = and £'(xh) > behaves as J 2 — > 0~. We have found that the 
function £ has no nodes, when U is determined by the fundamental nonabelian solution (for all 
allowed values of the parameters). Thus we conclude that the fundamental nonabelian solutions 
are stable. When U is determined by the n-th excitation the function £ seems to have exactly 
n nodes (we have checked this up to n = 2), hence the n-th excitation has n unstable modes. 



For uj 2 = and 



a 



the function £, after making n oscillations, tends to zero at infinity 



which signals the existence of static perturbations (bifurcation points). 

One can apply the same method to examine the stability of the RN solution. This was done 
by Lee et al. [9] in the full EYMH theory. They showed that the pulsation equation for the 
Higgs field (which in the case of RN solution decouples from the pulsation equation for the YM 
field) has no unstable modes. Thus the question of stability of the RN solution does not depend 
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on the parameter (3 (strength of the Higgs coupling) and reduces for every (3 to the Eq. (41), 
where now 

U-^^l. (43) 

For given xh > 1, the existence of negative eigenmodes in this potential depends on a. For 
a > (*max(%H) there are no negative eigenvalues, hence the RN solution is stable. For a < 
amax(jff)) the RN is unstable. Equivalently, one can say that for given a, the RN is stable if 
X H > x H ax (a), otherwise it is unstable. 

We have found numerically that when a decreases the RN solution picks up additional 
unstable modes. That is, there is a decreasing sequence {a n }, where ao = a mSLX , such that 
in the interval a n+ i < a < a n there are exactly n unstable modes (we have checked this up 
to n = 4). This is consistent with the fact that the RN solution has infinitely many unstable 
modes for a = 0, as follows easily from (43). We are convinced that the sequence {a n } coincides 
with the sequence {a^} of the bifurcation points at which the n-th excitation appears. Due 
to highly unstable behaviour of solutions near a™ ax we were able to verify this assertion with 
sufficient numerical accuracy only up to n = 1. However, as long as there are no other bifurcation 
points (as we believe), our conclusion is based on the general theorem of the bifurcation theory 
which says that if the operator governing small fluctuations if self-adjoint (hence its eigenvalues 
are real), then at the bifurcation point one eigenvalue passes through zero, whereas elsewhere 
the eigenvalues cannot change sign [10]. 

One could have anticipated the stability properties of solutions from a mere comparison of 
masses. In the region of the plane {o.,xh) where two distinct solutions coexist (with the same 
xh), the solution with lower mass (i.e. nonabelian) is stable, while the solution with higher mass 
(i.e. RN) is unstable. However, it should be emphasized that in more complicated situations 
the naive comparison of masses is not conclusive for stability. Actually, this happens in the full 
EYMH theory. To see this, let us consider the bifurcation diagram (m, a) for xh > 1 in the full 
EYMH theory. As follows from the results of Ref. 5 in the Prasad-Sommerfield limit and our 
preliminary results for finite f3, this diagram is more complicated than that shown in Fig. 3. We 
sketch it in Fig. 5. 

The horizontal line in Fig. 5 represents the mass of the RN solution. There are two branches 
of (fundamental) nonabelian solutions: the upper branch AB and the lower branch CB. These 
two branches merge at the bifurcation point B, which corresponds to a max (xn)- There is a 
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second bifurcation point A, corresponding to some ao(xn), where the upper branch merges 
with the RN solution. One can use these facts to infer the stability properties of solutions. As 
we have pointed out above the existence of a bifurcation point is a necessary condition for the 
transition between stability and instability. However, it is not a sufficient condition, since it is 
only when the lowest eigenvalue passes through zero that the stability changes. Therefore, if 
there are two branches which merge at the bifurcation point and one of them is stable then the 
other one will have exactly one unstable mode. Applying this reasoning in the present context, 
we see that if the lower branch is stable for small a (what we expect), then the whole lower 
branch must be stable, whereas the upper one is unstable with exactly one negative mode. By a 
similar argument the RN solution becomes unstable at «o- Therefore, in the region between ao 
and a max there are two stable solutions and one unstable. Note that in this region there is no 
relation between stability and mass. In particular the point in Fig. 5 at which the lower branch 
crosses the horizontal line is not related to change in stability (contrary to the suggestion in 
Ref. 5). 

The analogous bistability region exists in the plane (m, xh) for fixed a. This means that for 
given mass m there are two distinct stable black hole solutions which violates the weak no-hair 
conjecture [11]. 

6 Discussion 

We have analyzed static spherically black hole solutions in the strong Higgs coupling limit of 
the EYMH theory. We have found that the spectrum of solutions depends on the value of the 
dimensionless parameter a characterizing the model. Our results may be summarized as follows. 

For a > 1 there is only one branch of solutions, namely the RN family (parametrized by the 
radius of the horizon th > y/G/e). For a < 1 there appears a second (fundamental) branch of 
essentially nonabelian solutions. The radius of the horizon of these solutions is confined to the 
interval r™ n (a) < rn < r^ ax (a). The upper bound r^ ax (a) is a bifurcation point at which the 
two branches merge. The lower bound r^ m (a) is nonzero only for l/y/2 < a < 1. The solutions 
with rjj = r™ n (a) > are extremal in the sense that they have a degenerate event horizon 
and therefore their Hawking temperature vanishes. The analysis of linear radial perturbations 
about the solutions shows that the nonabelian solutions are stable whereas the RN solution 
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changes stability at the bifurcation point, i.e. it is stable for rjj > r^ ax (a) and unstable for 
r H < Vfj^^a). For sufficiently small values of a there exist also other branches of nonabelian 
solutions, which may be viewed as excitations of the fundamental solution. All excitations are 
unstable. 

These results have interesting implications for the fate of evaporating black holes in our 
model. Consider the RN black hole with unit magnetic charge and large mass. There are three 
different scenarios of its evolution due to the Hawking radiation depending on the value of a. 
For a > 1 the situation is the same as in the Einstein-Maxwell theory. The RN black hole emits 
thermal radiation and loses mass. When the horizon shrinks to rjj = VG/e (extremal limit), 
the temperature drops to zero and the evaporation stops. 

For a < 1, the RN black hole contracts to ru = rg ax (a) where it becomes classically 
unstable. The further evolution proceeds along the classically stable nonabelian branch and 
depends on whether a is smaller or greater than l/y/2. For a > l/y/2, the temperature drops 
to zero when the horizon shrinks to r™ n (a) (extremal limit) and the black hole settles down 
as an extremal nonabelian solution. For a < 1/V2, the temperature grows as the black hole 
contracts, so the black hole evaporates completely leaving behind a magnetic monopole remnant. 
This last scenario was first suggested by Lee et al. [6] as a possible consequence of the instability 
of RN solution. 
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Figure captions 

Fig.l The solution w(x) for xh = 2 and a = 0.01 (solid line), a = 0.2 (dashed line), and 
a = 0.288 (dotted line). 
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Fig.2 The bifurcation diagram in the (a, xh) plane. Excitations are not included. 

Fig.3 Mass m (in units \j\[Ge) of the fundamental nonabelian solution as a function of a for 
xh = 2. The mass of the RN solution is turn = 1.25. 

Fig.4 The n = 1 excitations for xh = 1- The dashed line represents the n = 1 colored black 
hole solution. 

Fig.5 The sketch of the bifurcation structure of black hole solutions with xr > 1 in full EYMH 
theory. Excitations are not included. 



